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Global existence of radial solutions for general semilinear 

hyperbolic systems in 3D 

Silu Yint Yi Zhou^ 


Abstract 

We study the well-posedness of radial solutions for general nonlinear hyperbolic systems 
in three dimensions. We give a proof of the global existence of radial solutions for general 
semilinear hyperbolic systems in 3D under null condition, with small scaling invariant 
^2,i(]]j3) -yYg obtain a bilinear estimate that is effective to the hyperbolic systems 

which do not have any time decay. It allows us to achieve the boundedness of the weighted 
BV norm of the radial solution. 


1 Introduction and Main Result 


We consider the following general first-order semilinear hyperbolic system 


du 

dt 




i=l 


du 

dx^ 


Q{u,u), 


( 1 ) 


where u = {pi,--- is the unknown function of {t,x^ ,x'^ ,x^), pi,--- , pi are 

scalar unknown functions valued in R, and vi, - ■ ■ , Vm are vector unknown functions valued in 
R^. Ai, A 2 , A^ are constant matrixes and Q(u, u) = m), Q’'(u, u))'^, Qp = (Qg • • • , Qf), 

QA = (Qi, • • • 5 Qm)^ with each = (Qp, is quadratic in u. 

System 0 is equipped with the following initial data, 

t = 0 : u = uo{x), X € R^. (2) 


In order to state our result precisely, we need to make some assumptions. The first assump¬ 
tion is the following. 

[HI]: Cauchy problem ([1]) Q is rotational invariant. 

We know that many physical systems in three dimensions like compressible Euler equations, 
relativistic compressible Euler equations and compressible MHD equations etc, are invariant 
under rotation of coordinates. The goal in this paper is to consider the global existence of 
radial solutions to Cauchy problem © which is invariant under rotation of coordinates. 
Consider the corresponding hyperbolic system in one dimension 

du . du 

a= (3) 

where u, Ai, Q{u, u) are defined as before in ([1]). 

We recall the following concept. 
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Definition 1.1. The system (0) satisfies null condition if any small plane wave solution u = 
u{t + is a constant in K.J to the linearized system 


du . du 

a +-"‘1? = “ 


is always a solution to the original system m- 


Then we assume 

[H2]: System (|31) satisfies null condition. 

In this paper, we only need to assume a weaker null condition for the corresponding one¬ 
dimensional system instead of the original system (HD- Roughly speaking, the assumption [H2] 
makes sure that waves of the same family do not interact in the nonlinear terms. 

The last but not least important assumption is presented as 
[H3]: The initial data satisfies 


ll'*^o|lvi/2,i(R3) < e, (4) 

where e is a small constant. 

Let us take a look at the invariance group of the system ©■ Suppose that u{t,x) satisfies 
system m, obviously u^{t,x) also satisfies system ([T|), which is defined by 

u^{t, x) = Xu{Xt, Ax), 


for all A > 0. That is, we have 


du^ . du^ 


dt 




i=l 


Differential (l5|) with respect to A and then take A = 1, we have 

i—1 

Here Lq denotes the scaling operator that introduced by Klainerman [lOj 

Lq = tdt + x^dj, 


( 5 ) 


and we use Einstein’s summation convention for repeated indices. Clearly, the norm given in 
0 is dimensionless that is scaling invariant 

II^o|Iu^2,1(r3) = 11^0 

As a primary task, we study on the coefficient matrixes under assumption [HI]. We have 
the following lemma, the proof will be given in section 2. 

Lemma 1.1. If system (Q]) satisfies [HI], then we system (Q]) should satisfy 
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and 

I m 

i,j—l s,q—l 

m 

Qliu^u) = '^fq^PjVq + Y X Vq, 

l<j<! S,o=l 

l<g<m 

fc = 1, • • • , 1 , p = 1, ■ ■ ■ ,m. 

The main result is given as follows. 

Theorem 1.1. Under assumptions [H1-H3], the Cauchy problem (QJ) (0) has a unique global 
radial solution 

p{t,x) = {pi{t,r),--- ,pi{t,r))'^, v{t,x) = {-wi{t,r),-■ ■ ,-w™(t,r))^, (6) 

r r 

where x = {x^,x'^,x^), r = 

We discover that some linear combinations of the unknown functions pi, Wi should satisfy 
wave equations, with nonlinear terms like rdr[(pa + Wa)iPb + wi,)], rdt[{pa + Wa){pb + Wb)], and 
{pa + u!a)ipb + yjb) where a,b € {1, • • ■ ,l+m}. Null condition ensures that a ^ bin the nonlinear 
terms. If we take notation F to represent these nonlinear terms, then by the following linear 
estimate that we will present in section 4 

Y [\\^o9tiru{t, •))||li(R+) + \\I^o9riru{t, •))||li(R+) 

a^0,l 

<c^holl 1T2,1 (r3) + E w^oFWl i([0,t]xR+); (7) 

a^0,l 
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one should give a suitable bilinear estimate on F to obtain the global existence. Once the bilinear 
estimate is proved, we complete the proof of Theorem 11.11 by the global iterative method. It is 
a standard process that can be followed by |22] . so we omit it here. Here and everywhere in 
this paper, notation C reprensents for some constant independent of the initial data. 

Theorem 1.2. Let i = a,b, consider the following wave equations 

I {dt - >Hdl){rpi) = Fi{t,r) 

]t = 0:rp^= rpio{r), {rpi)t = -Xi{rWio)rr{r), 


f (dt - Xjdl){rwi) = G^{t,r) 

[ t = 0 : rwi = rwioir), {rwi)t = -Xirpio{r), 

where wt is defined as 


Wi{r) = / Wi{s)ds. 

do 

If Xa Xb and XaXb 0, then there exists 


X! ll-^ 0 ^^r-[(Pa + Wa){pb + ^b)] HiHfO.t] xR+) 
a=0,l 

<C'^||pao|lvv2U(R3) + Il'^WaO II vi/ 2.1 (R3) + ^ ( II ^0 || Li ([O.t] xR+) + II ^0 || Li ([o,t] xR+)) ) 

a=0,l 

• (l|Pbo|lvv2,i(R3) + ||'Mibo|lvv2,i(R3) + ^2 (II^O.^b|lid[0,t]xR+) + ll^t'O l|Li([0.t] xR+)) ) 


• ( 8 ) 


a=0,l 


If one eigenvalue is equal to zero, without loss of generality, suppose Xa = 0, then pa, Wa 
satisfy ordinary differential equations 

{ Pat — fla 
^at — /*2a; 

and the bilinear estimate has a little difference. 


Theorem 1.3. Consider 

Pat — I^a (^T) 
t=Q ■. Pa= Paoir), 

{dt - XldDfrpb) = Fb{t,r) 
t = 0 : rpb = rpbofr), {rpb)t = -Xb{rwbo)rr{r), 
f {d? - Xldl){rwb) = Gb{t,r) 

[t = 0 : rwb = rWbo{r), {rwb)t = -Xbrpboir), 

then, we also have 

X! \\I^0'rdr[{Pa + 'Wa){Pb + yJb)]\\Lmo,t]xR+) 
a^0,l 

— C'^||Paollvi/2.i(R3) + Il'W'aO Ilvi/2,1(R3) + 'y ^ (l|-^'0^'5i'-^a||Li([0,t]xR+) + II ^0 II ([0,t] xR+) ) ^ 

O! = 0,l 

■ (l|Pbo|lvi/2.i(R3) + ll'*^bo|lvt/2.i(R3) + ^ (||Lo“^’b|Ui([0.b]xR,) + ||L?G,|Uq[0,b]xR,))). (9) 

a^0,l 


Wat = Ga{t,r) 
t = 0 : Wa = Wao{r), 
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Remark 1.1. Theorem, \1.2\ and Theorem fO are natural generalizations of the following well- 
known bilinear estimate in one spaee dimension 


“ 1 “ — -Ti(t^x') 

t = 0:(j)i = </>q(x), 


i = 1,2. // Ai ^ A 2 , then 

\\M2\\ Li([0,t]xR) < C'(ll<^olUi(R) + 11-^1 llii([0.t]xR)) (ll</'olUi(R) + ll-^2||Li([0.t]xR))- 
But the three dimensional estimate is much more difficult than the case of one space variable. 

It is worthy mention that the solution u{t, •) is not in the space for any time t > 0, 

although uo € rather in some weighted Sobolev space 


[W^odtiruft, •))|| li ( r +) + \\L^driru{t, •))|| li ( r +) 

a=0,l 


< Ce. 


Let us mention the existence theory of nonlinear wave equations given by Klainerman m- 
He proved the global existence of classical solutions for 3D scalar quasilinear wave equation 
based on two basic assumptions, smallness of the initial data and null condition of the non- 
linearities. This crucial work was also obtained independently by Christodoulou [4] using a 
conformal mapping method. However, both of their work depend on the decay of correspond¬ 
ing linear wave equations. In this paper, we develop this to the hyperbolic systems which do 
not have any time decay. 

Our long term project is trying to understand the formation of singularities and the global 
existence of weak solutions for systems of conservation laws in three dimensions. The multi¬ 
dimensional characteristics are so complicated that the tools which form the basis of a theory for 
hyperbolic conservation laws in a single space dimension do not extend to higher dimensions. 
Toward this goal, a key intermediate step is to understand solutions with radial symmetry 
exactly as that pointed out by Bressan et al in their book [3] . Unfortunately, radial symmetric 
systems do not admit global BV solutions generally (see O HH]), the theoretical analysis of 
several space dimensions still remains a grand challenge. Fortunately, we discover that the 
semilinear hyperbolic system allows global radial solutions with a bounded weighted BV norm. 
Thus to research the difficult hyperbolic systems in three dimensions, we consider the case with 
radial solutions as a first step. 

Another motivation for doing this paper comes from the one dimensional quasilinear hyper¬ 
bolic system 


Ut+ A{u)ua: =0, (10) 

where A(u) is an n-th order square matrix. Let Ai(u), • ■ • , A„(u) be the eigenvalues of A(u), and 
li(u), • • • , ln(u) denote its corresponding left eigenvectors. According to the important formula 
on the decomposition of waves which was given by John [^, \i Wi denotes the i-th component 
of Ux in (fT0|) , 


Wi = li{u)Ux, 


then we have 


Wit + {Kiu)Wi)x = '^lijk{u)WjWk. ( 11 ) 

Once we regard Xfu) and ^ijk{u) as constants, dill) is a semilinear system satisfying null 
condition. 


5 



For the completeness and for our future work, we also study the problem of quasilinear 
system in three dimensions. Consider the following general first-order quasilinear hyperbolic 
system 


du ^ A / \ n 


dt 


Z=1 


( 12 ) 


where u = (pi,-- - ,Vm)’^ is the unknown function of {t,x^,x'^,x^), pi,--- ,pi are 

scale unknown functions valued in R, and vi, - ■ ■ , Vm are vector unknown functions valued in 
Ai{u), A 2 {u), A^{u) are {I + 3m) x {I + 3m) matrixes with suitably smooth elements of u. 
We have 

Theorem 1.4. If the general three dimensional quasilinear system m is invariant under 
rotation of coordinates, and pk{t,x) = pk{t,r), Vp(t,x) = -Wp(t,r) are its radial solutions, then 
there exist bkq (p, w) and Cpj (p, w) such that 


Pkt + ^ hp{P: w)Wpr + ^ -bkp(p, w)Wp = 0 






(13) 


Wpt T ^ ^)Pjr 0- 

i=i 


Before ending the introduction, let us mention some important related works on classical 
solutions theory of quasilinear hyperbolic systems. In one spatial dimension case, the first result 
was given by John [6], he studied the formation of singularities when the system is genuinely 
nonlinear in the sense of Lax [12]. Liu [16] generalized John’s result in a neighborhood of u = 0, 
in which some linearly degenerate characteristic fields are allowed while other characteristics 
are genuinely nonlinear. Hormander [^ reproved the result given in [^ and obtained the sharp 
estimate for the life-span of solutions. By introducing the concept of weak linear degeneracy 
in m, the global existence and the life-span of solutions to the general hyperbolic systems 
for small initial data were completed. More results on this problem have been established for 
linearly degenerate characteristics or weakly linearly degenerate characteristics with different 
smallness assumptions on the initial data, such as [ 111151123 ]. However, the the corresponding 
theory for general one-order multi-dimensional hyperbolic systems is still open, because of the 
complexity of characteristics. We also refer a part of outstanding works for the existence theory 
of special hyperbolic systems in multiple dimensions. For 3D nonlinear wave equations, John [7] 
and Alinhac [1] proved that the solution of Cauchy problem with sufficiently small initial data 
but disobeying null condition will blow up at a finite time. John and Klainerman [S] showed the 
almost global existence theory for nonlinear scalar wave equation. Sideris and Tu m developed 
Klainerman’s generalized energy method in m to obtain global existence under null condition. 
For compressible Euler equations, Sideris [H] considered the formation of singularities for a 
polytropic ideal fluid in 3D, where the initial data are constant outside a compact set. and a 
complete analysis of the effect of damping on the regularity and large time behavior of smooth 
solutions was done in [20] . By constructing special explicit solutions with spherical symmetry, 
Li and Wang m investigated the blowup phenomena for compressible Euler equations. 

The remaining part of this paper is organized as follows. We devote section 2 to calculat¬ 
ing the conditions which coefficient matrixes should satisfy, when the system o is rotational 
invariant. And we give the proof of Theorem 11.41 We establish the null condition in section 3, 
and system © is simplified under conditions [HI] and [H2]. In section 4, we derive the linear 
estimate which is needed. In section 5, we introduce a core bilinear estimate with non zero 
eigenvalues. Then we perform the bilinear estimate with zero eigenvalues in section 6. 
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2 Rotational invariant system 


In this section, we study the general semilinear hyperbolic system © which is invariant under 
rotation of coordinates, and we will give the proof of Lemma If .11 in the following. 

In the case of rotational invariance, we shall treat scalar functions and vector functions 
separately. We set 


A,, = 


Ci jji I ’ 


where A’' = (aj^.), = (b®^), C* = (Cp^.), D’’ = (d^^), in which every element is a three 

dimensional row vector, is a three dimensional column vector, and is a 3 x 3 matrix, 
i = 1, 2,3; j, fc = 1, • • • ,1; p,q = 1, ■ ■ ■ , m. To be specific. 


Jil3\ 


^pq ^pq ^pq 

pq ^p<7 

J231 rl^32 J233 

^pq ^pq ^pq > 


^jq ’ ^pk \^pk'^ ^pk^ ^pk) 5 ^pq I ^ 

Suppose that the quadratic term of o has the following general formation; 

I m 

Qi{u,u) = Y. T%p,p, + ^ {T%T)l,T)l)p,v, + ^ u 




^ ^sq 
0^22 
^ ^sq 


L!i23 




^<3<l 

l<q<m 


s,q—l 


q/c31 oA;32 q/c33 
^^sq ^^sq ^^sq ‘ 


yq-! 


Q 7 (w,u)= ^ ^ 


i<j<i 

l<g<m 


™ /nffi npf^ 

+ E 

Vflffl ^pa33 




s,q—l 


where Tf^-, ff/, and fc = 1, • • ■ , p = 1, • ■ • , m; a, /3, 7 = 1, 2, 3. 

Proof of Lemma ll.il Set O = (oi, 02 , 03 )^ = (o^) € is any orthogonal matrix with 
detO = 1, and let y = Ox. Under rotation of coordinates, pk{t,x) and Vp{t,x) become 


Pk{t,Ox), 0'^Vp{t,0x), 


(14) 


k = !,■■■ ,l,p=l, 


, m. For simplicity, we introduce two notations, 

O^H4(oV,--- 

d{0 = d^^,+d^^2+d^^3, 

Then take (fTTl) into system ([T]), the left side becomes 

3 

[v(t, Oa;)] 


/ p{t,Ox) \ 

^YO^[vit,Ox)] 

^ ' 1—1 


A^ B^ 


d^. 


pit, Ox) 
0'^[v{t, Oa:)] 


Thus we have 


pit,y) 

vit,y) 


E 


E A^o,j X: \ 


j=i 


E o[C'' 


a=i 
3 


Va=i 


x; o[D3]o^ 


a=i 
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pit.y) 

yit,y) 


(15) 
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By assumption [HI], before and after coordinate rotation, the unknown functions should meet 
the same equation. So we obtain the following relationships, 




( ^3 \ _ O/pl ^3 

\^pk ’ ^pk 1 ^pk / \^pk •> ^pk ? ^pk ) ^ 


dpij — Odpm^ 


!(oi) 0 ^, 

Since O is arbitrary, there are 


dpq — Odpq 


( 02 ) 0 ’ 


dpq — Odpq 


( 03 ) 0 ^ 



which gives Ai, A 2 and H 3 in Lemma ll.il 

The nonlinear part of © can be treated in a similar way. We have algebraic conditions as 
follows, 


(T 
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^ ^ph2 ^ph3^ 
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V ^ jQ ’ ^ jQ ^ ^ .1 q) ~ ^ 
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= 0 , 




f^sq 
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which gives that 


fnPf^ ^p313X / 0 o' 

npf^ Off2 j^p323 ^ _^piii 0 0 

Voffi 0ff2 Off^y 0 0 0 


I m 

= Y1 

i,j—l s,q—l 


Qliu,u) = ^fq^PjVq + Y ^ ^ 9 ' 

1<3<1 s,q—l 

l<g<m 


This has completed the proof of Lemma 11.11 □ 

In order to simplify the notations, set bjq = b}P, Cpk = c^l, dpq = dpf. Take ([ 6 ]) into system 
HI), and using Lemma [1.11 we can easily get the following result by calculating directly. 

Theorem 2.1. If system (If) is rotational invariant and has radial solutions pk(t,x) = pk(t,r), 
Vp{t,x) = jWp{t,r), then it can be reduced to 


2 ^ ™ 

Pt + BWr + -Bw = Y^ ^ijPiPj + Y^ ^sqWsWq 

i,j—l s,q—l 

Wf CPr — ^ ^ jqPj^Q'i 

l<J<i 

l<q<m 


(16) 


where p = (pi, • • • , pi)'^, w = {wi, 


^sq — 


111 
sq 5 




‘ )'^m) ? B — {bjq)lxmj C — {ppk)mxl7 ^ij — 
qZ \T 'Y' _ /'Y'^11 '^mll\T A 'Y'm'v 

. i £J , L jq — y L ‘ ^ L J — [ 1 , 1 


^sqJ 


1 

ij^ ’ ^ ij' 
\T 


’ ^ 27 / 


Next, we will give a beginning of our future work and prove Theorem ll.di 
Proof of Theorem 11,41 Similarly, let 


Ai{u) = 


'A^{u) B\u) 
C\u) D\u) 


A^ = {a)^{u)), B^ = (b}fu)), C* = (c;,(w)),71* = (d;fu)), 


* = 1,2,3; j,fc = I,- - • ,/; = 1, 


, m, and 


(w) = {b]q{u),b)^{u),b]{{u)), cl^{u) = iBp^{u),c)p^{u),Bpi{u)) , 


= 


^dpMu) 4 f(ti) 

dfqlu) dfYiu) dffiu) 

di^u) df^iu) diHu) 


As the discussion in the semilinear case, if system (jl 2 [) is rotational invariant, then for any 
rotation matrix O, Ai{u) should satisfy 


a]kiP^v) 

a%{p,v) I = O 


^%{p,v) 


= O 


/^bl (p, 0 '^[u])^ 

k 2 / 


b^fP,O^H) 




i^pk{P 


a]fc(p,O^H)' 
a%ip,0^[v]) 

'-'jkiPA)/ \ajfc(P,O^H)/ VbJgiP,'-*" 

, Cpfe (p, u), (p, u)) = O(c^fe (p, [u]), (p, [u]), (p, [u]))O- 

dpf P, v) = 0(dif p, 0^[v])oii + d^f p, 0^[u])oi2 + d^f p, 0^[u])oi3)0' , 
dpfp,w) = 0(difp, 0^[v])o2i +d^fp, 0^[u])o22 Td^fp, 0^[u])o23)0^, 
dpfp,^^) = 0(difp, 0'^[u])o3i Td^fp, 0^[u])o32 Td^fp, 0'^[u])o33)0'^. 


(17) 
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Then we consider the radial solutions of the quasilinear system x) = pk{t, r), Vp{t, x) = 
:Wp{t, r), and take them into (TT^ . we have 


i 3 


m 3 




j = l 0!=1 
m 3 


q—1 a,0 — 1 

m 3 


X 
— 'I 

r 


Y. Y -bkqiP, y)Wq = 0 

q—1 a^0—l q—1 a—1 

13 m3 

Wpt + YY't^pj^p^^^p^-^Y Y ^ 2 '^pq y)y’qr 


j=l a=l 
m 3 


q—1 Oi,0—1 


-YY ^^^pqMWq+YY7^Pq^P^^'^^'i=^^ 


q—1 cx.,0—1 
m 3 


q—1 a=l 


where fc = 1 , • • • , = 1 , • • • , m, and = (c^^-, = {d^gP,d°^P, d^^^Y. 

Multiply (y)^ on both sides of the last equation, and take notice of (fTTl) . we have 

Y~^kjip^'^) = Y ^^^Kqip^y) = bll{p,o'^[v]), 


Q!=l 

3 


a,0—1 


J2 ^c'‘j(p,v) = ci](p,o^iv]), Y. = 


a,0—1 

So we can reduce (ITS)) to 


OL,0,y—l 


(18) 


m m 

PM + E hl\{p. 0^[v])wgr + E H\{p^ 0^[v])wg = 0 

9=1 9=1 

Wpt + E = 0. 

i=i 


(19) 


Because O is arbitrary, there exists bkq,Cpj such that bkqipjw) = 6 ^(p, 0^[u]),Cpj(p, w) = 
cl]{p,0^[v]). □ 


3 Null Condition 


We already know that Ai and Q{u, u) in ([3]) have special formats in Lemma I lTI In this section, 
we will study the nonlinear terms in (1161) under condition [H2]. If we rearrange the order of 
unknown functions in system ([3]), let V^°‘^ = {vf, • • • , v'^Y, a = 1,2, 3, then system ([3|) equals 
to the following 
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D 

0 

Ol 

1/(2) 

\v^^Y 


^0 

0 

0 

-d) 


i^i/(3) y 


/ QP{u,u) \ 
Q^'^iu, u) 
Q^'^ lu, u) 

\Q^^{u,u)J 


(20) 


where B = (&,,) e C = {cpk) e D = {dpq) e f)., = (()““, • • • 

and 


Qp = iQ{,---,Q?), Q™ = (Qr,--- ,Qr)- 
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Set 


Ml 


0 B 
C 0 



Let A = diag{Ai,--- ,Ai+m} be the eigenvalue matrix of Mi, L = (lij)(j_j^rn.)x(i+m) be its 
corresponding left eigenvector matrix, and R = (i?y )(;+m)x(i+m) = L~^. Similarly, let LM 2 R = 
A. These assumptions are reasonable because m is a hyperbolic system. 

For any fixed Ag, e € {1, • • • ,l m}, consider the plane wave solutions to the linearized 
system of ( 1 ^ . 


Z+m Z+m 

Pk — ^ ^ ^ ^ ^Z+p,a^^(^ 

a—1 a=l 

2m 2m 

Vp = ^ ^ Rpl)Yl)(^X Ae^), Vp — ^ ^ Rrn~\-p,bYl)(^X Ae^), 
b^l b^l 

in which 

Z m 

Ya Ag^) ^ ^ ^o^jPj ^ ^ q i 

q^l 

m m 

Yl)(^X Ag^) — ^ ^ Rq'^q “t“ ^ ^ R,m-\-q^q • 
q^l q^l 

Because of assumption [H2], we have 


z 

i,j=i 


m 


^ ^ ^sq^l+s,eRl+q,e — 0: 

(21) 

T^g-Rje-Rz+q,e = 0, 

l<j<l 

l<q<m 

(22) 


where k = 1, ■ ■ ■ ,1; p = 1, ■ ■ ■ ,m. 

We remark that if the system is not strictly hyperbolic, without loss of generality, suppose 
that Ai = A 2 , thus we furthermore obtain 


E ^^JR^lRJ 2 = E = 0 , 


i,j=l 




(23) 


for all fc = 1, • • • ,1. 

Next we pay attention to (HU). Rewrite L = G ^ii+m)xi^ pi ^ ]^(i+m)xm^ 

multiply block on both sides of the first part of (HU) and multiply block on both sides of 
the second. Then we have 

2 ^ z+m Z m 

Pet ~t“ AgUZgT- -|- “AgUZo, — ^ ^ lek ^ ^ ^ ^ ^ijRiaRjb “t' ^ ^ ^sq^i-\-s,aRl-\-q,b 

k—1 a,b—l i,j—l s,q—l 

< ■ (Pa + Wa){pb + Wb) 

m Z+m 

'^et H“ ^ePer — ^ ^ R,l-\-p E E ^ jqRjaRl-\-q,b{Pa H“ '^a)(Pfc H“ '^b) ■, 

p—1 a,b—l i<i<i 

I, l<g<m 
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where 


i 

Pa — ^ ^ ^aj Pj: '^a 


m 

la,l+qWq. 

9=1 


Take notice of (ED) and (ESI, we have 


Pet “t“ AeXCer “t“ Agdie — ^le 


ihei “t“ ^ePer — /*2e 


(24) 


(25) 


in which /le and / 2 e are linear combinations of {pa + Wa){pb + Wb) and e = 1, •••,/ + m. It 
is crucial that a ^ hhy the algebraic condition (ED and (1221) . Furthermore, if the system has 
repeated eigenvalues Ai = A 2 , then 5^2 (a ^ 2) when a = 1 (6 = 1 ) because of (l23)l . 


4 Linear Estimate 

In the following contents of this article, we always use an uniform notation /e to represent for 
/le and / 2 e appeared in (1^ . We have 


{d1 - Xld^){rpe) = rdtfe + rdrfe + fe 




T 

'){rWe) ='^ j 


dtfedr + r fe¬ 


rn 


But if Ae = 0, we should better replace (E^ by 


Pet — fie 
Wet = f2e- 


(27) 


Take notation F as the nonlinear part of the first equation of (051) . If Ae 7 ^ 0, we know that 

df - Xldl = {dt ± Xedr){dt T Xedr), 


and 


[Lo, DAe] = -2nAe, 

in which Da,, denotes the wave operator with wave speed Ae- Then 

{dt + Xedr){dt - Xedr)L^{rpe) = (Tq + 2)“F, 

{dt - Xedr){dt + Xedr)L^{rpe) = (Tq + 2)“F, 

where Lq = tdt + rdr, a = 0,1. The right side of the above equation can be written as LqF 
briefly, this will have no essential difference in our analysis. 

Let 


A = {dt- Xedr)Lo{rpe), B = {dt + Xedr)Lo{rpe), 

then we have 

{dt + Xedr)\A\ = sgn{A)L'^F, 

{dt - Xedr)\B\ = sgn{B)LQF. 
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Thus 


II^I|l1(R+) < II^o||l 1(R+) + ll■^O^IUT[ 0 .^]xR+) < C'I|wo|IvV2.i(B 3) + II^O^IUi([ 0 .i]xR+), 
II^IIl1(R+) < II^o||l 1(R+) + 11-^0-^llLi([ 0 ,t]xR+) < C'l|wo|lvh2.i(R3) + II^O-P"llLi([ 0 ,t]xR+)- 


That is 


||9tL^(rpe) + 9rTo(’'Pe)IUi(R+) ^ C'll^io|lu/2.i(R3) + ||io^|lLi([0,t]xR+)- 






If Ae = 0, it is much easier to get (l2^ because of and [Lo,dt] = —dt- We do not repeat 
that process here. 

In order to obtain the global existence, by contraction mapping principle and global iterative 
method, we hope to have a uniform bound for pe and Wg and their derivation with some weighted 
norm. The key to this is the bound for the quadratic terms. 

5 Bilinear Estimate with non zero eigenvalue 

Even if the nonlinear terms of (1261) have forms like rdtfe, rdrfe and fe, we only need to estimate 
rdrfe- Spirits of the proofs are similar, but the estimate of rdrfe in some suitable weighted 
Sobolev space is much more complicated than others. 

Proof of Theorem 11.21 We would estimate the bilinear estimate under the case of Aa = 
—Xb separately. First, we pay attention to |Aa| 7 ^ |Ab|. 

For simplicity, according to Duhamel principle, we can only consider the following equations 



where i = a,b. We would only study for a = 0 because of [To,n] = —20. Thus, we need to 
perform the proof of the following 


Ikar[(pa + Wa)(Ph + W{,)]||il([o,t]xR+) <C'(^||s(pa||Li(R+) + || (s'l/'a)'|1 Li (R+)) 

• (l|s^(¥’h)'||Li(R+) + ||s(sV'b)"llLbR+)) ■ (29) 


Since the wave speeds are positive and equals to |Ai|, we can always assume A^ > 0 for notation 
convenience. If (pi and tpi are spherical symmetric, we can express pi and Wi in the form (see 


i) 



then 



(30) 
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where sh^{s) = 2 ^[s(pi(s) + (s' 0 i(s))'], rni{s) = ^' 0 i(s)- 
Notice that 


rdr[{pa +Wa){pb + W&)](t,r) =dr[r{pa + Wo)] ' {pb + Wb) - 2{pa + Wo) ' [Pb + Wb) 

+ (Pa + Wa) ■ dr[r{pb + Wfc)] 

=Ia+Ib+Ic- (31) 


We first estimate I a- Let Eq be small enough, we reduce matters to considering two cases: 
Case 1: r > Sot. 

Case 2: r < Eot. 

Our strategy for dealing with 1/r is as follows. There is nonsingularity when r is away from 
the origin, so we can treat is directly. If r is near to the origin, we only need to contribute one 
derivative by Hadamard’s formula. 

Case 1: By (1301) we have 



\lA\drdT < 


{XaT + r)ha{\aT + t) 


(•XbT+r 


- / shb{s)ds 

^ J\XbT-r\ 

\bT+r 


+ [ [ \Kt - r\ha{\XaT - r\)- f shi,{s)d. 

do oM_|_ ^ d\XbT—r\ 


'0 ^R+ 
rt 


1 


f AaT+r 


|A|,T-r| 

Ai,r+r 


, {sma{s)yds 

^ d\X.r-r\ 


+ 


1 


r* AaT+r 


'|Aar-r| 


i{s)ds 


J|AfcT-r| 

XbT+r 


shb{s)ds 


' |A 6 r-T-| 

fXbT+r 


shb{s)ds 


n^bT+r 

{XaT + r)ha{XaT + r)— / smb{s)ds 

;_l_ ^ J\XhT—r\ 

pXbr-\-r 

IAoT - rl/iodAoT - r|)— / sm;,(s)d, 

:+ ^ •llAi.r —r| 


drdr 

drdr 

drdr 

drdr 

drdr 


+ 


+ 


iQ ./R+ 
r-t 


f'XaT-\-r 


dAbr-rl 
f AfcT+r 


^ A|A„r-r| 


{sma{s)yds / smb{s)ds 


1 


r* AaT+r 


'|A6T-r| 

fXbT+r 


drdr 

drdr 


^ A|A„r-r| 


sma{s)ds 


MAtr-rl 


smb{s)ds 


drdr 


—Iai + Ia2 + I A3 + Ia4 + Ia5 + I AO + Ia7 + Iao- 


Let us considering Ia 2 and Iao and then other terms are easy to follow, we estimate it in 
two subcases. For both cases we shall need to set ^ = XaT — r. 

Subcase 1: epr < r < e(d^|CI- Under this condition, we have ^ t < £9 

from the integral expression of Iao, Ave know that 0 < s < (-^ + l)r. Hence, 


Ia2 <||s^h||Li(R+) 


|(AaT - r)ha{\XaT - r|)| ■ -drdT 


0 •teoT<r<Eo ICI 


-dr 

(|y+l))-b«l ^ 


<U||shdUi(R+) [ l^haiOm 

«/M-l_ 

f* f \ 1 

/a 6 <||smb||L=o(R^) / / \{XaT - r)ha{\XaT - r\)— / ds\drdT 

♦A 0 ^ SoT<r<€Q^\^\ ^ ^0 
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-dr 


f /-eo^lfl 

<C'||(sm,)'|Ul(R^) / |C/ia(5)Me / 

Jm+ J(^+i)-i|4l r 

<C'||(sm 6 )'||Li(R^)||C/la||Li(B+)- 

Subcase 2: r > Eot, and r > £(7^l^l- case, 

|Aa - Abir < AblAaT - r| + AalA^T - r| < XbSor + \a\XbT - r|. 

Take notice of Aa 7 ^ Ab, AaAb 7 ^ 0, so Eq is sufficiently small and should be less than 
have 

< |Abr - r| < s < AbT + r < (— + l)r, 

^Aa eo 


2\i 


, we 


which shows 


£0 ^ ^ 2Ao 

s < r < — -:—rS. 


Ab + £0 


lAa-AbT 


Hence, 


Ia 2< ^ \{XaT-r)ha{\XaT-r\)\-\- ° shb{s)ds\drdT 


r>eg |« 

2 Aa 

<c [ / \shb{s)\ds [ *’ -dr 

Jr+ Jr+ ^ 


^b+gp , 


2 Aa 


2 Xa 

'l>a->bl® 1 

-IT dr 

=o_„ 


<C'||^^a||Li(R+) ||s/ib||Li(R+)) 

Ia6< f I ^ |(AaT - r)/ia(|AaT - r|)| ■ 1^ [ “ smb{s)ds\drdT 

<c [ \^haiO\d^ [ \smbis)\ds 

<C [ \^ha{0\d^ [ \smbis)-\ds 

Jr+ Jr+ s 

<C'IIC^a||Li(R+)||TOb||Li(R^). 

If we estimate Iai and Ia 5 when r > Eqt, we should set rj = XaT + r. There naturally 
holds r < £q ^ry as long as £0 is small enough. We omit the specific calculation here. And if we 
integral I as, Iaa and I at, I as by parts with respect to r, then it is same as Iai, Ia 2 , I as or 
I AO- There is no essential difference with them. 

Case 2: Let ^0 = XaT + 9r where 6 G [—1,1], since r < Eqt, then 

^ 6 < r < 


Ao + £0 


Xa £0 


We rewrite the integral of \Ia\ as 


|/.4|drdr < 

/O 7R+ Jo JR+ 

rt 


{XaT + r)ha{XaT + r) ■ {pb + drWb) 


drdT 


0 JR+ 


{XaT - r)ha{XaT - r) ■ {pb + drWb) 


drdT 
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1 


^AaT+r 


^ IAtiT—r| 

^ /*Aar+r 


(sTOa(s))'ds • (pb + drWb) 


J\\aT-r\ 

=Ia9 + ^^10 + Iaii + Ia12- 


sma{s)ds ■ {pb + drWb) 


drdr 

drdr 


Using — L°“ estimate, which can be traced back to the 1960s, see Lemma 1 in reference m, 
we have 


Thus, 


|pb + drWb\ < Ct (||r Pb|Ui(R+) + ||r V’"IIli(R+))' 


<C'(||r2(PbllLi(R+) + IIli(b+)) [ [ \{XaT+ r)haiXaT+ r)-\drdT 

Jo Jr<eoT ^ 

1 

<C{\\r^^b\\LHm+) + IkVb IIli(R+)) / l6^a(Ci)M6 / ° ° -dr 
<C||6^a||Li(R+)(|k^<PbllLi(R+) + Ik^V'b IIli(R+)), 


Iaii <C{\KPb\\mR+) + Ik V'b IIli(k+)) 


-|[(Aar + r)ma{XaT + r) - {XaT - r)ma{XaT - r)]i|drdT 


<C'(lkVbl|Li(K+) + IkVh IUi(K+)) / / / 

Jo Jr<£oT J — 1 

<C'(lk^<^hllLi(K+) + IkVh IIli(k+)) / de [ |(sma)'| , MCe 


dO—drdr 

T 


s—Xar-\-9r 

1 


—dr 


Aa+eo 


^9 


<C\\{i9maY\\L^R+)i\\r^‘Pb\\L^R+) + IkVb IIl1(R+))- 


Obviously, the estimates of Iaio, Iai 2 are also easy to follow because of Ia 9 , Iaii and the 
property of integration by parts. The method of treating Is and Ic can be imitated as Ia- 
There is no essential difficulty, so we omit the detailed proof here. 

Had not forgotten that, we have one last mission in this section, i.e. to prove the estimate 
when Xa = —Xb = X. Without loss of generality, we assume that A > 0. According to the above 
proof, we discover that only subcase 2 in case 1 is highly dependent on |Aa| 7 ^ |Ab|. Thus, we 
should separately consider the case for r > Sqt and r > where = Xt — r. Fortunately, 

we can get more information from the solutions of the following linear equations 


r {df - X^dfKrpa) = 0 

[ t = 0 : rpo = rpao, dtrpa = -X{rWao)r 

( (52 - X^d^,){rpb) = 0 

11 = 0 : rpb = rpbo, dtrpb = X{rwbo)rr 

then we can prove that 


/ id! - X^d!)ir^,) = 0 
[ t = 0 : rwa = rwao, dtrwa = -Xrpao, 

{d! - X^d!){rwb) = 0 

t = 0 : rwb = rwbo, dtrwb = -Xrpbo, 


||r5,.[(p„ + zi}a)(p6 + Wh)]||Li([0.t]xR+) <C'(lkkaolUhR+) + II ™a0 lU^ («+) ) 

■ (l|Pbo|lvv2,i(R3) + llwbollu/^’HRS)^ 


(32) 
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Using D’Alembert formula, we have 


r{pa + Wa) ={r - Xt)pao{r - Xt) + Wao{r - Xt) + (r - Xt)wao{r - Xt) 

1 ~ ~ 1 
-—[{r - Xt)wao{r - Xt) + {r + Xt)wao{r + Xt)] + — spaois)ds, (33) 

J I At—r| 


r{pb + Wb) =(r + Xt)pbo(r + Xt) + Wbo{r + Xt) + (r + Xt)wbo{r + Xt) 

1 ~ - 1 
- + Xt) + {r - Xt)wbo{r - Xt)] + — spbo{s)ds. (34) 

2?' 2r Jixt-r\ 

Then 

dr[r{pa + Wa)] ={r - Xt)p'^g{r - Xt) + pao{r - Xt) + 2wao{r - Xt) + (r - Xt)w'^g{r - Xt) 

- 7 ^[(r - Xt)wao{r - Xt) + Wao{r - Xt) + (r + Xt)wao{r + Xt) + Wao{r + Xt) 
Zr 

+ (r - Xt)pao{r - Xt) - (r + Xt)pao{r + At)] 

^ j'Xt+r 

+-^[{r - Xt)wao{r - Xt) + [r + Xt)wao{r + Xt) - j s/5ao(s)(is]. (35) 

J\\t-r\ 

We just show the proof of following three terms of I a = drlr^pa+Wa)] ■ (P6 + W{,), and the others 
can be estimated similarly. Since r > Eqt and r > we have r < p = r + Xr < (l + A£(C^)r. 

Inserting (IMl) and (1551) into I a, each term can be followed as one of the following 

Pao{\r - Ar|) • Wbo{]r - Ar|)-^^^^— Zzl^drdr 

I 

<C [ / |Pao(|^|)who(|CI)4'M’’'^^ 

v/b+ Jr>£g bsl ^ 

<C'||who||L»(R+) [ |Pao(|CI)l|c^’’C^^ 

Jr+ s 

<(^11^001^1(8 +) I|W{,oIIli(R+)i 



r/. 

r>e 

-Wao(k - 

eo-r 1 j. 

Ar|) • -Wboii 
r 

+ Xt) I dr dr 


r>eybei 

r — At|) Idrdr'j ( f f 

^ ^Jo Jr>{l+ 


(/j^' Wao(s)ds)2 


(/o’'w6o(s)ds)2 

V 


/* 1 ^ 1 ^ p p'n ^ 

Kois)dsdyyy wlo{s)dsdpy 

<(^(^ l?|wao(?7)t^C) " Vwlo(p)dsdpy 


<(^ll(CWao)'||Li(R+)||(?yWbo)'||Li(R+), 


|—w^Q(r + Ar)|drdT 


1 

2 



nXt-\-r ^ nXt-\-r 

/ (xpaoia)da— j 

'|Ai —r| ^ olAi—r| 


spbo{s)ds\drdT 
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drclrdsda 


<c / \apao{(T)\ / |spbo(s)| 


<C'l|Pao||Li(R+)||P6o||Li(R+)- 


r>(i+^)-io 


''r<en 


We deal Ib and Ic with the same method like I a when Aq = — Ah = A. By now the proof 
of Theorem ll.2l is finished, since ||s/ih||Li(R+) < ll'S^/ihl|Li(R+)■ □ 


6 Bilinear Estimate with zero eigenvalue 

In this section, we consider situation with zero eigenvalue. 

Proof of Theorem 11.31 By Duhamel principle and along the discussion in section 5, we 
need to estimate 


|k9^[(Po + Wa)(p& + Wt,)]||Li([0.t]xR+) 

<C'(lkPaollihK+) + ll™lollLhR+)) ■ (ll'5^‘/^blUhK+) + l|s('S'0b)”llLi(R+)) ■ (36) 

From (1551) . we know that Aq and Ah play different roles. This makes us think of dealing 
with Aq = 0 and Ah = 0 differently. Certainly, it is impossible that both of them equal to zero 
because of the null condition (| 2 T]) and (l 2 ^ . 

If Aa = 0, Ah > 0, then pa, Wa satisfy ordinary differential equations, 


Pat — fla 
djat — f2a- 


(37) 


Using expression OT, we have 

nXbr+r 

\r{p'ao+dao)- shb{s)ds\drdT 

(+ ^ ^IAfcT—r| 

/•AbT+r 

\r{Pao + Wao)— / smb{s)ds\drdT 

(+ ^ t/IAbT —r| 


—^Ah “t” IAtu- 


They can be estimated as 

Iah < f [ \r{p'^Q + w'^Q)-[HbiXbT + r) - Hb{XbT - r)]\drdT 

Jo Js.+ ’’’ 

- [ [ \'>’ip'aO + d'ao) f (XbT + 9r)hb{XbT + 6r)d9\drdT 

Jo Jr+ J -1 

<C ( \r{Pao + d'ao)\dr j \£,hb{(,)\dr, 


I Am ^ / \^{Pa 

Jr+ 

<ess sup 
rGR+ 


I 


aO w Wao 


p2 


XbT+r 


J\\bT-r\ 

<C'lk(p;,o+«^lo)llii(R+) , 


p\bT+r 

' -^|A|,T-r| 

smb{s)dsdT 
1 , 


K+ Jo 


smb{s)dsdT\dr 

Hp'ao + dao)\dr 
drdr 


„ , smb(s)ds 

J\Xbr-r\ 
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lk(PaO + Wao)IUi(R+) 


rt I 1 r^bT+r 


(smb(s)yds drdr 


r^hT+r 

3 / {smb{s)yds drdr. 

1 + ' J\XhT—r\ 


dAml — 


‘ /■ I 1 


0 ^R+ ' ' J\\tT-r\ 


{smb{s)yds drdr 


= I I \\{XbT + r)mb{XbT + r) - \XbT - r\mb{\XbT - r\)\drdT, 

Jo Jr+ ^ 

we estimate it by considering two cases: 

Case 1: r > Eot. In this case, we have ^ = A^r + r < + l)r, 77 = A^r — r < XbT + r < 

(|r + I)?-, then 


f f ^\{XbT + r)mbiXbT + r)\drdT <C f f ^l^rubiOldrd^ 

Jo Jr>eoT ^ Jr+ d)'>(^ + l)-i{ ^ 

<C [ \mbiO\d^, 


[ [ \\XbT - r\\mb{\XbT - r\\)\drdT <C [ [ ■\\r]mbir])\drdr] 

Jo Jr>€Qr ^ Jr+ J+ ^ 


<C [ \mb{r])\dr]. 

JRa 


Case 2: r < Eqt. Then = XgT + 9r > XgT — r > — l)r, and 

lAmi < / {smb{s)y'\s=Xi,r+0rd9 drdr 

«/0 Jr<eQT J ~1 

<C r f ( \{^gmb{£,g))"\drd^ed9 <C f M^gmbm)"\d^0- 

J -1 Jr+ dR+ 


ll-fB||Li([o,t]xR+) < / / \{PaO+u!ao)- shb{s)ds\drdT 

Jo Jw^ ^ J\XbT—r\ 

/-AfoT+r 

\iPaO+Wao)— smbis)ds\drdT, 

’+ ^ J I A^r —r I 

then the estimate of ||/B||Li([o,i]xR+) is like ||/yi||Li([o.*]><»+) = 0, Ab 7^ 0. 

We calculate dr{r{pb + Wb)), 

dr{r{pb + Wb)) =[(A6T + r)hb{XbT + r) - {XbT - r)hb{XbT - r)] 

2^ 2^ r^bT+r 

-[(AfcT + r)m6(A{,T + r) - (AfcT - r)m{,(A6r - r)] + — / smb{s)ds, 

^ J\XbT-r\ 

thus we can estimate ||dc||Li([o,t]xR+) by the method what we have used as above. 
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If Aa > 0, Ab = 0, then pb, Wb satisfy ordinary differential equations, and pb{t,r) = pbo{r), 
Wb{t,r) = Wbo{r). From the previous discussion, we only need to pay attention to the case of 
f < Sq'T- We know 

rdr-iiPa + Wa){pb + Wb)] =dr[r{pa + Wa)] ’ {Pb + Wb) + {pa + Wa) ■ rdr{pb + Wb) 

- {Pa+Wa){pb+Wb), (38) 


and expression 


2^ pXbt+r 

Pb + Wb=- / 

^ J\\bt-r\ 


sgb{s)ds + dr J smb{s) 


s—\bt-\-9r 


dO 


—Xbt-\-9r 


9d9, 


1 f 2 ) 

where pb = --^Pb i will help us to discuss. 
We obtain 


n {XaT ± r)ha{XaT ± r) ■ {pbo{r) + Wbo{r)) 

<SoT 

<C ( [ \^±ha{^±) ■ {pbo{r) + Wbo{r))\drd^± 

<C'll^±^o(^±)||Li(R+)||pbO + W6 o||l1 (R+), 

[ f f g^9d9-{pbo{r)+Wbo{r)) 

Jo Jr^soT J— 1 ^ 

f i f li^emai^e))'■ {pbo{r) + WbQ{r))\drd^ed9 

J-1 Jr+ 

<C'||(^ema)'||Li(R+)||pbO + Wbo||Li(R+), 


drdr 


drdr 


where = XaT ztr, = XaT + 9r. Thus 

rt 


<eo^ 
t r /.I 


drdr 


n {Pa + Wa) ■ rdr{pb + Wb) 

<eoT 

^9a{s)\^^^^^^g^d9 + J ^{smay\^^^^^^g^9d9 ■ r{p'go + w'bo) 


<C 

/O JM.+ J-1 

<C'[||Ce5a||Li(R+) + ||(Cem-a)'||Li(R+)]||rp(,o + ™bollii(R+) 


(39) 


drdr 


<C 


/ / {Pa+ Wa) ■ {pb + Wb) 

Jo J r<eoT 


<eo^ 
t r /•! 


drdr 


d0 + y (sma)'l 


U=A r+9r^^^ ' iPbO + Wbo) 

/o -'IR+ J -1 “ ■ J -1 

<C'[||Ce5a||Li(R+) + ||(CeWa)'||Li(R+)]||pbO + Wbo||Li(R+)- 


drdr 


And other terms are easier to perform. We have completed the proof of Theorem 11.31 □ 

Acknowledgement The authors would like to thank Prof. Gui-Qiang Chen for the en¬ 
lightening discussions and encouragements. 


20 














References 


[1] Alinhac, S., Rank 2 singular solutions for quasilinear wave equations. Internal. Math. Res. 
Notices, 18(2000), 955-984. 

[2] Bressan, A., Contractive metrics for nonlinear hyperbolic systems. Indiana Univ. Math. 
J., 37(1988), 409-421. 

[3] Bressan, A., Chen, G. Q., Lewicka M., Wang, D. H., Nonlinear Conservation Laws and 
Applications. The IMA Volumes in Mathematics and its Applications 153, Springer, 2011. 

[4] Christodoulou, D., Global solutions of nonlinear hyperbolic equations for small data. 
Comm. Pure Appl. Math., 39(1986), 267-282. 

[5] Hormander L., The lifespan of classical solutions of nonlinear hyperbolic equations. Lecture 
Notes in Math., Springer, 1256(1987), 214-280. 

[6] John, F., Formation of singularities in one-dimesional nonlinear wave propagation. Comm. 
Pure Appl. Math., 27(1974), 377-405. 

[7] John, F., Blow-up for quasilinear wave equations in three space dimensions. Comm. Pure 
Appl. Math., 34(1981), no. 1, 2951. 

[8] John, F., Klainerman, S., Almost global existence to nonlinear wave equations in three 
space dimensions. Comm. Pure Appl. Math., 37(1984), no. 4, 443-455. 

[9] Klainerman, S., Solutions to quasilinear wave equations in three space dimensions. Comm. 
Pure Appl. Math., 36(1983), 325-344. 

[10] Klainerman, S., Uniform decay estimates and the Lorentz invariance of the classical wave 
equation. Comm. Pure Appl. Math., 38(1985), no. 3, 321-332. 

[11] Klainerman, S., The null condition and global existence to nonlinear wave equations. Lect. 
in Appl. Math., 23(1986), 293-326. 

[12] Lax, P. D., Hyperbolic systems of conservation laws 11. Comm. Pure Appl. Math., 10(1957), 
537-567. 

[13] Li, T., Wang, D. H., Blowup phenomena of solutions to the Euler equations for compressible 
fluid flow. J. Dijf. Eqs., 221(2006), no. 1, 91-101. 

[14] Li, T. T., Zhou, Y., Kong, D. X., Weak linear degeneracy and global classical solutions for 
general quasilinear hyperbolic systems. Comm. Partial Diff. Eqs., 19(1994), 1263-1317. 

[15] Li, T. T., Zhou, Y., Kong, D. X., Global classical solutions for general quasilinear hy¬ 
perbolic systems with decay initial data. Nonlinear Anal., Theory, Meth Appl, 28(1997), 
1299-1322. 

[16] Liu, T. P., Development of singularities in the nonlinear waves for quasilinear hyperbolic 
partial differential equations. J. Diff. Eqs., 33(1979), 92-111. 

[17] Marshall, B., Strass, W., Wainger S., — L‘^ estimates for the Klein-Gordon equation. J. 

Math. Pures Appl., 59(1980), no. 9, 417-440. 

[18] Rauch, J., BV estimates fail for most quasilinear hyperbolic systems in dimensions greater 
than one. Comm. Math. Phys., 106(1986), 481-484. 

[19] Sideris, T. G., Formation of singularities in three-dimensional compressible fluids. Comm. 
Math. Phys., 101(1985), no. 4, 475-485. 


21 



[20] Sideris, T. C., Thomases, B., Wang, D. H., Long time behavior of solutions to the 3D 
compressible Euler equations with damping. Comm. Partial Diff. Eqs., 28(2003), no. 3-4, 
795-816. 

[21] Sideris, T. C., Tu, S. Y., Global existence for systems of nonlinear wave equations in 3D 
with multiple speeds. SIAM J. Math. Anal., 33(2001), no. 2, 477-488. 

[22] Sogge, C. D., Lectures on Non-linear Wave Equations. International Press of Boston, Inc., 
2008. 

[23] Zhou, Y., Global classical solutions to quasilinear hyperbolic systems with weak linear 
degeneracy. Chin. Ann. Math. Ser. B, 25(2004), 37-56. 


22 



